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A sheet of elastic foil rolled into a cylinder and deformed between two parallel plates acts as a non- 
Hookean spring if deformed normally to the axis. For large deformations the elastic force shows an 
interesting inverse squares dependence on the interplate distance [Siber and Buljan, arXiv: 1007.4699 
(2010)]. The phenomenon has been used as a basis for an experimental problem at the 41st Interna- 
tional Physics Olympiad. We show that the corresponding variational problem for the equilibrium 
energy of the deformed cylinder is equivalent to a minimum action description of a simple gravi- 
tational pendulum with an amplitude of 90°. We use this analogy to show that the power-law of 
the force is exact for distances less than a critical value. An analytical solution for the elastic force 
is found and confirmed by measurements over a range of deformations covering both linear and 
non-Hookean behavior. 



I. INTRODUCTION 

In a recent study Siber and Buljan analyze the follow- 
ing simple yet pedagogically rich problem from the theory 
of elasticity. A thin flat elastic sheet (e.g., a piece of plas- 
tic foil) is rolled into a cylinder (radius bo) and placed be- 
tween two impenetrable plates which are parallel to each 
other and to the axis of the cylinder, see Fig. [TJ The dis- 
tance 2b between the plates is fixed externally. For b < bo 
the foil acts as a spring exerting a force of magnitude 
F(b) on each of the plates. An interesting property for 
this kind of spring is the non-Hookean power-law scaling 
of the elastic forced F oc b~ 2 , which holds (as we show 
here, exactly) for b < b c ~ 0.7&o- Measuring F(b) in this 
universal scaling regime has been proposed 1 as a method 
to determine bending rigidity for such objects as plastic 
foils, electrical connectors, biological membranes and mi- 
crotubules, possibly nanotubes and monolayer materials 
(e.g., graphene) etc. A lab problem based on measuring 
of F(b) for standard plastic transparency films has been 
recently given to world's top secondary school physics 
students at the 41st International Physics Olympiad (Za- 
greb, Croatia, 2010). 

The corresponding mathematical problem of con- 
strained minimization of the foil's elastic energy may 
seem difficult and hardly illuminating. Existing solution 1 
is based on analytical approximations and numerical fi- 
nite element optimization while a textbook approach- 
relies on force equilibrium conditions for strongly bent 
elastic rods that are rarely covered in standard physics 
curricula. 

In this paper we formulate and solve the variational 
problem using tangential angle parametrization of the 
profile shape. This reveals equivalence to another con- 
ceptually rich physics system — the large-amplitude pen- 
dulum: elastic energy of the foil maps onto kinetic en- 
ergy of the pendulum while the fixed inter-plate distance 
translates into the cosine-shape potential energy. Simple 
mechanical considerations allow us to deduce the exact 
inverse squares law for the elastic force at b < b c . Using 
the standard solution for the large-amplitude pendulum 




FIG. 1. A thin-walled elastic tube deformed between two 
parallel plates. The thick black contour line marks the profile 
at the critical value of b c /&o = 0.717770, thinner lines above 
and below correspond to different b/bo between and 1. 

in terms of elliptic functions, we find (a) the exact values 
for b c and related constants; (b) a single transcendental 
equation that determines F(b) for b > b c ; and (c) com- 
pact analytic form of the universal profile. Finally, the 
deduced functional dependence F(b) is compared to mea- 
surements on a plastic film in the table-top setup used 
in Ref. 1. This show feasibility of a quantitative demon- 
stration of both the universal non-Hookean law for b < b c 
and the usual linear regime for b — > bo. 

II. FORMULATION OF THE PROBLEM 

Specific property of the deformation geometry in this 
problem is "cylindricity" : the Gaussian curvature van- 
ishes at every point. This eliminates non-uniform stretch- 
ing/compression (§14 of Ref. d) and leaves only the bend- 
ing contribution to the total elastic energy W e \ for thin 
sheets (thickness d <C bo). The problem is essentially one- 
dimensional and the energy functional ([1]) is the same as 
for an elastic filament (Kirchhoff rod, see §18 of Ref. 0): 

W el = ^ [lC 2 ds. (1) 
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FIG. 2. (a) Coordinate system for profile energy calculation, 
(b) Equivalent pendulum problem with 90° maximal devia- 
tion, (c) A possible four contact-point profile. 



(This effects of gravity are ignored.) Here k = (1 — 
v 2 )~ 1 Ed 3 /12 is the bending rigidity, E is the bulk Young 
modulus of the foil material, v is the Poisson ratio, h is 
the length of the cylinder in the non-deformed direction 
(i.e. parallel to the axis), and K. is the planar curvature 
of the deformation profile in the plane normal to the axis 
(x-y plane). The shape of the profile is described us- 
ing natural parametrization^ {x(s),y(s)} in the coordi- 
nate system defined in Fig. [2£i with origin at point O; 
s is the arc length measured counter-clockwise form O. 
Integration in ([1]) is performed along the entire profile, 
and the absence of stretching^ implies that J ds = 2nbo 
regardless of b. The constraint imposed by the plates is 
expressed by requesting that the vertical coordinate of 
point O' is 2b (see Fig.^i), 



x(s = irba) = 2b . 



(2) 



In the following we find the minimum of W e \ subject to 
the constraint and taking into account impenetrabil- 
ity of the plates. The elastic force is then obtained from 
F(b) = -dW e i/d(2b). 



III. ANALYTIC SOLUTION 

A. Lagrangian formulation and analogy to the 
pendulum 

An intrinsic quantity characterizing the shape of a pla- 
nar curve is its tangential angle 0(s) defined as {i, y} — 
{cos 0, sin 0} (dot denotes derivative with respect to s). 
Particular advantage of employing 0(s) for the present 
problem is that the (signed) curvature equals^ to simply 
JC = 8. Integrating x = cos 8 over the right half of the 
profile between the contact points O and O' gives the 
expression of the constraint ^ in integral form: 



l cos0(s) ds = 2b = — — / ds 
Jo ™o J 



(3) 



At the points of contact between the foil and the plates, 
the tangent to the profile must be parallel to the y-axis, 



hence the boundary conditions for 0(s), 

0(s = O) = -7T/2 , 0(s = 7T& O ) = +7T/2 . (4) 

A standard way of turning a constrained optimization 
problem to an unconstrained one is the use of Lagrange 
multipliers. In our problem both the target function W e \ 
and the constraint ([3]) are expressed as integrals (func- 
tionals) of the unknown function 8(s). Denoting the La- 
grange multiplier for Eq. ([3]) by Wq (for reasons that will 
become clear shortly) and using left-right symmetry in 
Eq. ([T]), the variational problem becomes that of uncon- 
strained minimization with respect to a function 8(s) and 
a number u> 2 of the following functional 



2 Kh 



tWj 



c 



0(s),<?(s),w o 2 



ds 



where 



C = -9 2 + uj 2 cos 8 - w, 



2b 



(5) 



(6) 



Now the analogy to the pendulum problem has become 
manifest. Viewing W as dynamical action and s as time, 
one recognizes in Eq. ([5]) the Lagrange function of a rigid 
pendulum in uniform gravitational field with the angular 
frequency of small oscillations equal to loq (see schematic 
drawing in Fig. Eb)- In the pendulum problem, 8 is the 
angle of deviation from the stable equilibrium which must 
satisfy the Newton's equation (Euler-Lagrange equation 
of the variational problem)^, 



uj n sin f 



0. 



(7) 



In standard mechanics problems ujo is usually given, and 
the fulfillment of the boundary conditions Q is ensured 
by choosing the appropriate initial velocity 0(0) = /Co 
which is a function of ujq. In our case, both ujq and /Co are 
not known a priori and must be determined by satisfying 
both the boundary conditions ([U and the constraint (J3J) . 
The optimal value of luq is proportional the elastic force 
itself since, after "integrating out" the pendulum degree 
of freedom [i.e. substituting 0(s) in the Lagrange function 
with the solution to the equation of motion (|7J|] , W(ujo, b) 
satisfies dW/duj — and thus 



F 



-nh 



dW{uj 0l b) 
db 



= ~ K h~Qj^~ = 2k/icjq . (8) 



We shall use this neat property of Lagrange multipliers 
for the calculation of F(b). 

Before tackling the problem of finding /Co and loq 
(which is mostly mathematical), we develop physical ar- 
guments for the universal scaling behavior which consti- 
tute, in author's opinion, the most important message of 
this article. 
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B. Universal regime, b < b c 

Qualitatively, /Co (6) starts from the value of /Cq(6o) = 
1/&0 and decreases as b is decreased. At a ceratin crit- 
ical (in the sense of separating two qualitatively differ- 
ent behaviors) b = b c the contact curvature vanishes, 
K-o(^c) = 0- The corresponding pendulum problem be- 
comes that of free oscillations with amplitude 7r/2 and 
frequency ujo(b c ) = u> c ("critical pendulum"). Period of 
these oscillations is well-known^ [also derived below as 
a special case of Eq. (flB])]. T = 4K(l/2)/u c , where K 
is complete elliptic integral of the first kind. Since ac- 
cording to Eq. (Q} it must be equal to To = 2/T&o, we 
get 

w c = (Co bo)' 1 with Co = K -i= ' = 0.847213 .... (9) 

v 7r 

When b is decreased further below b c , an extra condi- 
tion not accounted for by the Lagrangian formulation ([5]) 
becomes relevant: the plates do not allow the foil to bend 
outwards, thus /Co remains zero also for b < b c . In order 
to accommodate the imposed small values of b without 
violating the rigid plates, a finite part of the foil in the 
vicinity of s = and s = irbo must remain flat, /C(s) = 
(these parts are marked by horizontal dashed lines in the 
profiles shown in Fig. [T]). Concurrently, the sections of 
the profile that do bend obey Newton's equation ([7]), al- 
though with some u>o(b) > uj c . Assuming continuity in 
IC(s), we conclude that the deformed part of the profile 
for b < b c must start with /Co = (same as for b = b c ) 
but cover a length shorter than the full length 2nbo of 
the foil cross-section. This shortening of the deformed 
part can only be accommodated by a faster pendulum 
since the frequency uj a (b) of the latter remains the only 
adjustable (i.e. 6-dependent) parameter in the problem 
for b < b c . (The other parameter, /Co, is pinned to zero 
by the presence of a finite flat part.) But, according to 
Eq. ([7]), changing wo results in a uniform rescaling of 
the arc length parameter s, therefore the bent parts of 
the profile at b < b c (marked by curved lines between the 
plates in Fig.[l} must be geometrically similar to the cor- 
responding halves of the critical profile at b = b c (marked 
with a thick contour in Fig. [T]). 

Having established that the critical profile shape is also 
universal (in the sense of being applicable to any b < b c ) 1 
we can determine the scaling law for the force using sim- 
ple dimensional considerations. Scaling of b and bo with 
b/bo fixed does not change the overall shape of the pro- 
file — the tangential angle 9 remains the same function 
of s/bo- Therefore the elastic energy ([1} can be written 
a&i W c \(b) = nhU(b/bo)/bo where U is a dimensionless 
function of b/bo- For b < b c the flat parts of the pro- 
file do not contribute to the energy. The curved parts 
are similar to the critical profile but b/b c times smaller, 
thus their energy must have a fixed U{b/bo) — > U(b c /bo) 
and rescaled bo — > bob/b c . This argument proves exact 
scaling for the energy, W e \ cx b , and, consequently, the 
force, F cx b~ 2 for b < b c . The prefactor in the scaling 



law can be described as follows: 

F(b) = F c § = 2^M^Y 2 for b < b c . (10) 
So 

where F c = F(b = b c ) and Eqs. © and © have 
been used. The phenomenological "stadium profile" 
model, considered as a variational ansatz in Ref. U, as- 
sumes circle as a universal profile and predicts F(b) = 
(■jr/2)nh/b~ 2 , similar to Eq. (fT0|) . The exact critical 
shape is, however, more efficient than a circle in mini- 
mizing the bending energy. We can calculate the numer- 
ical prefactor in Eq. (p~0|) exactly by determining b c . To 
this end, we proceed to the integration of the pendulum's 
equation of motion. 

The first integral of Newton's equation is the energy 
conservation law, 9/2 — Wgcosf? = const. Taking into 
account cos#(s = 0) = gives 

\o 2 = ^ICo 2 +uj 2 q cos 9. (11) 

Integrating both sides of Eq. (fTTj) with respect to s and 
with respect to 9, then using Eq. ([3]) gives 

W= fj 9 2 ds = ^ + 2u 2 b, (12) 
o 

_ 1 r+*/ 2 

W=- 9d9 = ojoIo(1Co/u Jo ), (13) 

1 J-tt/2 

where 

, + 77/2 

I n (a) = (1/2) (a 2 +2 cos 9) 1 ' 2 - n d9. (14) 

J-tt/2 

Now b c can be determined from Eqs. (|T2"1) and (| 13)) : 
substituting b c , uj c , and for b, u>o, and /Co respectively, 
and using an identity Co = ^o(0)/2, one gets 

6 C = / (0)/(2w c ) = Co bo = 0.717770 . . . x b . (15) 

The results (fTCTI) and (|T5|) confirm phenomenological 
Eq. (17) of Ref. [l] (with numerical factor 0.912 there cor- 
rected to 4C^/?r = 0.91389). 

C. Small deformations, bo>b>b c 

In this regime /Co > and the relation between the 
period and the boundary conditions becomes more com- 
plicated. Expressing loq ds via d9 from Eq. (jlip and in- 
tegrating over from s = Q to s = irbo, we obtain the 
required additional equation, 

Trb oJo = 2Ii(Kq/uiq) . (16) 

This equation will provide the answer for the force, 
thanks to Eq. ([H]). The remaining unknown is the pa- 
rameter a = /Cq/wq, for which a single transcendental 
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FIG. 3. (a) Graphical solution of Eq. (JTTJl. (b) The curvature 
at the contact point /Co as a function of b/bo. The dashed 
straight line (b — 6 c )/(6o — 6 C ) is plotted for visual guidance. 



equation is easily derived by combining Eqs. (TT2"j) . (TT5T) . 
and fTS]). 
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7T / I (a) 
4 



(17) 



The functions /„, denned as integrals and appearing 
in Eqs. (fl"3j) . (|T6|) and (fl7l) . can be reduced to stan- 
dard elliptic integrals^ Io{a) — (4/fc)£(7r/4, k 2 ) and 
2i(a) = fcJ-"(7r/4, fc 2 ), where £ and J 7 are incomplete 
elliptic integrals of the first and the second kind, re- 
spectively, and k = 2/y2 + a 2 is the elliptic modulus. 
7i(0) = K(l/2) = tt/(2Co) in consistence with Eq. 
We were unable, however, to solve the transcendental 
Eq. (|17|) analytically and present a graphical solution in- 
stead, see Fig. [3^. 

As 6/60 changes from b c /bo to 1, the only root of 
Eq. (fT?) goes from to 00. Divergence of a — K-o/wo 
as b — > bo is consistent with ui 2 oc F 0. The curvature 
at the contact point /Co = au>o is plotted separately in 
Fig. [5b; /Co decreases from l/b to as 6 goes from bo 
to 6 C , in accordance with our previous qualitative discus- 
sion. 

Summing up the solution for b > b c , the exact elastic 
energy and the force are 



F = 



4nh 

—Io{a)h{a) , 
2kH (2I x (a) x 1 



(18) 
(19) 



where a(b/b ) is the root of Eq. (fTT)) . see Fig. [5^. 



D. Stability of the universal solution and 
bending-induced tension 

After formulating an explicit solution for b > 6 C , 
Eq. (|18l) . we can justify the universal scaling for b < b c 
more rigorously. It has been suggested^ that at b < 0.3&o 



profiles with a slight curvature in the horizontal part may 
provide a better solution than a "stadium" with com- 
pletely flat segments. Such a four-contact-point alterna- 
tive solution is sketched in Fig. (the concave segments 
BC and AD are exaggerated). The following argument 
shows, however, that any four-contact-point profile will 
have larger bending energy than the optimal profile with 
straight BC and AD. 

Let us fix the positions of the contact points on the 
profile, e.g., by gluing the foil to the plates at A, B, C, 
and D, The curvature JCa must be the same at all of 
these four points due to symmetry; furthermore, K-a is 
evidently positive. Thus the strongly curved parts AB 
and CD must conform to our solution for /Co = JCa > 
with a smaller bo, adjusted to their respective shorter 
lengths. Consider now increasing slightly the lengths of 
AB and CD while making BC and AD shorter so that 
the total circumference of the profile remains the same. 
In this new configuration the energy of the segments BC 
and AD will evidently be reduced by shortening. How- 
ever, the corresponding lengthening of the segments AB 
and CD shall also reduce their energy: for fixed b, longer 
(higher bo) two-point profiles have lower minimal total 
energy. This can be shown explicitly using Eqs. ©, (|16l) . 
(EE! and (fill): 



Fn 



dW r] 



d(2irbo) 



K,h 



ICi<0. 



(20) 



b— const 



Thus one can always lower the bending energy of a four- 
point profile like the one shown in Fig. [2fc down to the 
point where BC and AD become completely flat and 
can no longer be shortened. This minimal energy limit 
is precisely our zero-contact-curvature universal stadium 
profile. 

It is instructive to note that Fx/h is, by definition, 
the surface tension energy of the deformed foil. Uni- 
form tension may be induced in a cylindrically deformed 
sheet by bending even if the stretching deformation is 
negligible^ The value (|2"0)) for Ft and its independence 
of s can be confirmed using a formal method of local La- 
grange multiplies- which takes local tension into account 
explicitly. 

Recalling that our sheet is assumed to be made of uni- 
form elastic material with bulk Young modulus E, we 
can now estimate the neglected stretching deformation. 
Bending-induced negative surface tension is balanced by 
uniform linear stretching. By using Hooke's law, we can 
express the increment in the profile circumference 27rA6o 
due to stretching as A60 = &o-Pr / '{Edh). The corre- 
sponding contribution of stretching A Wgtrctch ~ Ehd b /6g 
to the total energy W ~ Ehd 3 /b is negligible if d 2 <C b^, 
thus confirming our approximation of an unstretchable 
foil. The universal profile is tension free, Fx(b < b c ) = 0. 
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FIG. 4. Elastic force F as a function of the half-height b, 
scaled by the critical force F c and the undeformed cylinder 
radius bo, in linear and inverse square root representations. 
Circles: experimental data scaled by a single fitting param- 
eter, F c . The spring was first gradually loaded (o) up to 
frmin = 0.4 bo and then unloaded (•). 



E. Critical shape 

Before discussing our final results for the force, 
Eqs. (fTU)) and (JTHJ), we briefly comment on the shape of 
the critical profile. Using the explicit solution of the pen- 
dulum problem [obtained, e.g., by integrating Eq. (jTTJ) , 
see Ref.0 for a pedagogical presentation] , and using prop- 
erties of the elliptic functions^ the shape of the profile 
at b — b c , in the units corresponding to b = l/(£o\/2)) 
is given parametrically for 9 g [— 7r/2,7r/2] by 



x{8) =Co/V2- 
y{9) =±v / cos^ 



(21) 



with the arc length s{9) = tt6 /2 + F(6/2,\/2). This 
shape is shown in Fig. [T] by the thick contour touching 
the plates. Profile shapes for b > b c are obtained by 
integrating Eq. (J7J numerically with the initial condition 
obtained from Eqs. (fTB")) and (JTTJ) and shown by thin lines 
reaching above the upper plate in Fig. [T] 



IV. DISCUSSION AND COMPARISON TO 
EXPERIMENT 

The final results for the force are shown in Fig. @] We 
use F c = 2k/i/((o&o) 2 to scale the force into a dimension- 
less form, F(b)/F c , which is a single universal function 
of b/bo- This function is shown by continuous lines in 
Fig.Hboth in terms oiF(b)/F c and (F{b)/F c )- 1 ' 2 to re- 
veal the region of power-law scaling. The point at b = b c 
is an inflection point of F(b) (with a jump in the third 
derivative) . 

It is instructive to verify that the rolled foil behaves 
as a Hookean spring when close to cylindrical shape (i.e. 
as b —> bo). Using the definition of I n (a) to get the 
large a expansion, we obtain a from Eq. (|17[) . bo — b ^ 
(tt/4 - 2/ir)a- 2 bo, and F(b) from Eq. (Q]JJ} as follows 



F(b) 



2nh bo - b An 



b 



for b 



■■>o ■ 



(22) 



This linear behavior is marked by a dashed line in lower 
right corner of Fig. 0J 

Measurements of F(b) have been performed using the 
setup and one of the samples (blue plastic binding covers, 
Set 1) described in Ref. [l|. The results are shown in Fig.0] 
by circles on a linear scale (cf. the logarithmic scale used 
in Fig. 5 of Ref. [fl). The spring was gradually loaded 
from b — bo down to & min = OAbo (open circles) and then 
unloaded back by increasing b up to zero force (filled 
circles). Note the hysteresis due to inelastic deforma- 
tions. The corresponding energy losses (hysteresis loop 
area) are 4% of W e i(& m in)- The critical force F c = 116 N 
with relative error of 3% (estimated from the residuals 
between the open data points and the analytical fit) cor- 
responds to the bending rigidity of k = (1.48 ± 0.04) mJ 
in reasonable agreement with Ref. [l|. 
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